USAAVLABS  TECHNICAL  REPORT  66-62 


to 


MICROMECHANICAL  BEHAVIOR 
OF  FIBER  REINFORCED  PLASTICS 


By 


Juan  Haener 


C  L  E 

FOR  FEDE 

TECH:'' 

\  R  1  N  G  H  0 

RAL  SCIEN1 
CAL  INFOR1 

Ha: dcopy 

%  i  cc 

Microfiche 

$  ,£  *> 

i  ft 

/ 

mw  (gf 

||  c 

September  1966 


‘n',^77 

\v 


FORT  EUSTIS,  VIRGINIA 


CONTRACT  DA  44-177-AMC-320(T) 
WHITTAKER  CORPORATION 


Distribution  of  this 
document  is  unlimited 


SAN  DIEGO,  CALIFORNIA 


D 


D  C 


fna?rarPDn  nrprp 
UJ 


DEC  2  2  1966 

Lb  Lj  Lb  Lj  L/  Lb] 

C 


OJUftOWJL  UAY  -  -  ; 


-h 

vaiT#. 

AL/QlLUSTaj* 


Di  sclaim  er  s 


The  findings  in  this  report  are  not  to  be  construed  as  an  official  Depart¬ 
ment  of  the  Army  position  unless  so  designated  by  other  authorized 
documents. 

When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a  definitely  related  Government 
procurement  operation,  the  United  States  Government  thereby  incurs  no 
responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way  supplied  the  said 
drawings,  specifications,  or  other  data  is  not  to  be  regarded  by  impli¬ 
cation  or  otherwise  as  in  any  manner  licensing  the  holder  or  any  other 
person  or  corporation,  or  conveying  any  rights  or  permission,  to  manu¬ 
facture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be 
related  the reto. 

Trade  names  cited  in  this  report  do  not  constitute  an  official  endorse¬ 
ment  or  approval  of  the  use  of  such  commercial  hardware  or  software. 


Disposition  Instructions 

* 

Destroy  this  report  when  no  longer  needed.  Do  not  return  it  to 
origin  ~tor . 


SIM  Nr 

»»ITE  SECT* 

I,  ^  WTF  SECTION 

mmmaf) 

WWCATIO 

. ./-tL..,  Qts* 

^wniin/imiiiuTr  coats  | 

••ST.  |  MIL.  rn/ar  iftCiA 


f 


DEPARTMENT  OF  THE  ARMY 

U.  S.  ARMY  AVIATION  MATERIEL  LABORATORIES 
FORT  EUSTIS.  VIRGINIA  23604 


This  program  was  carried  out  under  Contract  DA  44-177-AMC-320(T) 
with  Whittaker  Corporation,  Narmco  Research  and  Development 
Division. 


The  data  contained  in  this  report  are  the  result  of  research  con¬ 
ducted  to  determine  the  mechanism  of  load  transfer  through  a  com¬ 
posite  material  consisting  of  fibers  embedded  in  a  matrix  under 
consideration  of  Internal  stresses  caused  by  polymerization 
and  temperature  shrinkage  of  the  matrix  as  well  as  externally 
Imposed  loads. 

The  report  has  been  reviewed  by  the  U.S.  Army  Aviation  Materiel 
Laboratories  and  is  considered  to  be  technically  sound.  It  is 
published  for  the  exchange  of  Information  and  the  stimulation  of 


future  research. 


Project  1P121401A14176 
Contract  DA  44-177-AMC-320(T) 
USAAVLABS  Technical  Report  66-62 

September  1966 


MICROMECHANICAL  BEHAVIOR 
OF  FIBER  REINFORCED  PLASTICS 

Final  Report 

by 

Juan  Haener 


Prepared  by 
WHITTAKER  CORPORATION 
Narmco  Research  &  Development  Division 

San  Diego,  California 
for 

U.S.  ARMY  AVIATION  MATERIEL  LABORATORIES 
FORT  EUSTIS,  VIRGINIA 


Distribution  of  this 
document  is  unlimited 


OJUQOUL  HAY 


*■  I  i*  .  ^  4 


.  .'.XU  £UT± 

ADQXJAA 


_  IV 


ABSTRACT 


V 

The  distribution  of  stresses  in  s  unidirectional ly  oriented,  multi¬ 
fiber  composite  as  a  result  of  process  shrinkage  and  external  loads  has 
been  analysed  after  assuming  continuity  boundary  conditions  at  the  inter¬ 
face  and  certain  hexagonal  boundary  conditions  in  the  space  between  the 
reinforcements.  Selected  model  specimens  have  been  numerically  analyzed 
by  computer,  and  the  results  have  been  compared  with  photoelastic  experi¬ 


mental  studies 
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FOREWORD 


This  report  was  prepared  by  Whittaker  Corporation,  Narmco  Research  & 
Development  Division,  San  Diego,  California,  under  Contract  DA  44-177-AMC- 
320(T),  entitled  "Micromechanical  Behavior  of  Fiber  Reinforced  Plastics." 
The  work  was  accomplished  under  the  technical  supervision  of  Dr.  R.  Echols, 
Chief,  Physical  Sciences  Laboratories  Division,  US  Army  Aviation  Materiel 
Laboratories  (USAAVLABS),  Fort  Eustis,  Virginia. 

This  report  covers  work  conducted  from  14  June  1965  to  14  May  1966. 

The  principal  technical  investigator  on  this  project  was  Dr.  Juan 
Haener.  Dr.  Gerhard  Nowak  contributed  as  a  consultant,  and  Dr.  George* 
Burgin  was  responsible  for  computer  work.  Other  contributors  to  this  pro¬ 
gram  were  Messrs.  Krishna  Naik,  Ming-Yuan  Feng,  Noel  Ashbaugh,  and  Curt 
Thompson.  The  program  was  conducted  by  Narmco* s  Engineering  Department 
under  the  supervision  of  Mr.  B.  L.  Duft,  Engineering  Department  Manager. 

The  program  management  responsibility  resided  with  Mr.  B.  Levenetz, 
Research  Engineering,  and  the  administrative  responsibility  with  Mr.  R. 
Hidde,  Project  Engineering. 
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SUMMARY 


This  program  was  initiated  as  a  continuation  of  analytical  research 
on  composite  materials  which  had  been  reported  in  USAAVLABS  Technical 
Report  65-58.  The  objective  is  to  improve  the  understanding  of  the  mech¬ 
anism  of  load  transfer  through  the  composite  material  consisting  of  fi¬ 
bers  embedded  in  a  matrix  under  consideration  of  internal  stresses  caused 
by  polymerization  and  temperature  shrinkage  of  the  matrix  as  well  as  ex¬ 
ternally  imposed  loads.  It  is  anticipated  that  this  research  will  lead 
to  more  reliable  methods  of  predicting  mechanical  behavior  of  fiber  rein¬ 
forced  plastic  composites  under  load  in  order  that  the  structural  integ¬ 
rity  of  airframe  components  from  this  material  can  be  improved. 

During  the  program  reported  here,  the  boundary  conditions  at  the  in¬ 
terface  and  at  the  hexagon  boundary  have  been  analytically  investigated 
and  a  tridimensional  stress  and  strain  analysis  on  a  unidirectlonally 
oriented  multifiber  bundle  has  been  performed.  The  developed  analytical 
expressions  have  been  checked  on  a  special  case  of  a  single  fiber  and  then 
applied  to  analyze  the  distribution  of  stresses  within  a  multifiber  com¬ 
posite.  The  solutions  were  checked  by  computer  analysis,  which  also 
provided  numerical  results  for  stress  distribution  in  selected  model  spec¬ 
imens.  The  analytical  work  has  been  compared  with  results  of  microphoto¬ 
elastic  studies  conducted  on  similar  experimental  model  specimens. 
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INTRODUCTION 


The  high  strength-to-weight  ratio  of  oriented-filament  composites  can 
be  most  fully  realized  only  if  the  optimization  of  these  materials  is 
based  on  accurate  stress  analysis.  Test  results  have  shown  that  present 
stress  equations  do  not  accurately  reflect  the  true  stress-strain  distri¬ 
bution  in  matrices  reinforced  with  oriented  filaments.  At  present,  there 
is  incomplete  understanding  of  the  three-dimensional  mechanical  load 
transfer  between  reinforcements  and  matrix  in  a  composite  subjected  to 
externally  imposed  stresses  combined  with  internal  residual  stresses. 

The  influence  of  these  combined  stress  conditions  on  composite  character¬ 
istics  must  be  resolved  by  these  fundamental  steps: 

1.  Study  of  the  solutions  to  the  differential  equations  of 
distortions 

2.  Expression  of  the  physical  problem  in  terms  of  restricted 
numbers  of  these  solutions 

3.  Examination  of  mathematical  implication  of  the  solutions 

4.  Comparison  of  the  solutions  with  test  results 

5.  Derivation  of  practical  design  formulas  from  the  theoreti¬ 
cal  solutions 

During  the  process  of  the  previous  contract  [DA  44-177-AMC-208(T)] , 
mathematical  relationships  were  derived  for  a  single  fiber  embedded  in  a 
resin  cylinder  for  the  case  of  shrinkage  and  static  loading  with  general 
boundary  conditions.  Efforts  on  the  multifiber  case  have  led  to  the  point 
of  obtaining  solutions  to  the  differential  equations;  however,  boundary 
conditions  which  would  relegate  these  solutions  to  a  specific  physical 
problem  were  not  yet  applied. 

During  the  present  contract,  the  general  boundary  conditions  have  been 
specialized  to  correspond  with  the  experimental  specimens.  This  provided 
the  possibility  of  programming  the  equations  for  a  computer  analysis  and 
of  obtaining  numerical  results  for  stress  distribution  within  the  compos¬ 
ite  models. 

The  question  of  so-called  "potentials"  present  in  the  anaysis  of  the 
multifiber  case  has  also  been  resolved.  The  solutions  of  the  short  and 
very  long  multifiber  case  have  been  obtained  in  three  dimensions.  They 
have  been  programmed  for  the  very  long  fiber  bundle. 

As  a  special  case  of  the  multifiber  problem,  the  stress  distributions 
for  the  singie-fiber  problem  have  been  established  for  any  fiber  length. 
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Because  the  theoretical  analysis  is  applicable  for  any  combination  of 
elastic  properties  of  the  constituent  materials,  specific  properties  which 
correspond  to  those  of  the  selected  specimen  materials  have  been  used  in 
the  computer  analysis  to  establish  the  diagrams  for  stress  distributions. 
This  now  makes  it  possible  to  relate  the  theoretical  results  with  results 
of  the  microphotoelastic  studies. 


3 


TECHNICAL  DISCUSSION 


METHOD  OF  APPROACH 

There  are  two  generally  known  methods  for  reducing  the  fifteen  linear 
field  equations.  One  method  leads  to  three  equations  for  the  displacements 


dXj  G 


o. 

G 


*1 

at2 


(i) 


i  »  1,2,3 

j  -  1,2,3 


and  the  other  leads  to  six  equations  for  stresses. 

In  the  present  work,  the  displacement  method  expressed  in  equation  (1) 
was  utilized.  This  equation  covers  completely  the  linear  three-dimensional 
elastomechanics,  Including  the  mechanics  of  elastic  wave  propagations  (Ref¬ 
erences  1  through  5). 

One  of  the  reasons  for  this  decision  was  the  existence  of  extensive 
literature  (chapters  on  basic  equations  in  References  1  through  5)  showing 
how  facts  must  be  arranged  to  obtain  solutions  in  terms  of  a  restricted 
number  of  elementary  solutions.  Further,  it  was  believed  that  the  obtained 
solutions  of  the  displacement  equations  for  the  static  case  can  later  be 
more  readily  extended  for  dynamic  considerations  such  as  acoustical  and 
shock-wave  propagation  (Reference  4)  than  it  would  be  by  utilization  of 
stress  equations. 

Another  reason  for  selecting  the  approach  of  displacement  functions  was 
that  in  three-dimensional  elasticity,  this  approach  is  the  only  one  that  has 
been  used  effectively  in  obtaining  solutions  of  practical  problems. 

For  the  static  case  and  by  neglecting  the  body  forces  and  assuming 
rectangular  coordinates,  the  equation  for  displacements  becomes 


T3  V8  fi  -  0  . 


(2) 
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There  are  two  classical  approaches  often  used  (see  Reference  2)  to 
find  expressions  for  displacement  and  stresses  frpm  equation  (2):  the 
Galerkin-Westergaard  approach  and  that  of  Boussinesg-Fapkovitch  and  Neuber. 
Here,  the  second  method  was  selected,  because  the  stresses  appear  as  second 
and  not  as  third  derivatives  of  a  potential  P  defined  as  Papkovitch  func¬ 
tions.  This  fact  facilitates  considerably  the  transformation  into  cylindri¬ 
cal  coordinates  as  well  as  the  nunerical  calculations. 

The  displacements  expressed  in  Papkovitch  functions  are 

*i  "  ?i  '  4(1  -  v)  (Xj  PJ  +  *o)  (3) 

i  -  1,2,3 

J  "  1,2,3  . 


The  four  Papkovitch  potentials  in  equation  (3)  are  solutions  of  the  harmonic 
equations. 


V3  Pi  -  0  i  -  0, 1,2,3  . 


(4) 


The  completeness  of  the  Papkovitch  functions  is  shown  in  Reference  (6)  and 
also  in  Reference  (2) . 

Transformation  of  the  equations  into  cylindrical  coordinates,  a  process 
which  is  shown  in  Appendix  A,  equations  (Al)  through  (A25),  enables  us  to 
express  the  internal  elastic  condition  in  the  composite.  Without  going  into 
further  details  of  transforming  the  basic  equations,  (1)  and  (2),  and  equa¬ 
tion  (4),  we  transform  directly  the  solutions  and  the  F^  into  cylin¬ 

drical  coordinates,  obtaining 
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sin(n0)  +  Cn 
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(5) 


J  -  1, 2,3.0 


Equation  (5)  represents  three  Papkovitch  functions  for  each  of  the  two  mate¬ 
rials  with  36  constants  in  total.  Introducing  equation  (5)  into  equation 
(3)  and  after  transforming  (3),  we  obtain  for  the  displacements  (see  equa¬ 
tions  (A60)  through  (A76)  in  Appendix  A): 


in  radial  direction, 

§i  “  pi 


_ I _ 

4(1  -  v) 


»  (r  p‘ +  2  p° +  p°] » 


(6) 


in  tangential  direction. 


*•  ■  p»  -*o-h>  (7) 

in  axial  direction, 

5s  "  P3  ’  4(1  l  v)  ^  [r  P1  +  2  P3  +  Po]  •  (8) 


The  stress  distribution  from  strain  components  is  given  by  the  physical 
components  of  the  stress  tensor. 
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Here  Djj  are  the  elements  of  the  strain  tensor  in  cylindrical  coordinates 
and  are  restated  as  follows: 
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The  stresses  can  be  represented  and  calculated  directly  from  the  Papkovitch 
potentials  through  the  following  relation: 
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5ijV  9xt  9Xj  J 

*k  Pk  +  P°  j 

(16) 


i  -  1,2,3  j 


1,2,3  k  -  1,2,3. 


Equations  (5)  through  (16)  are  general,  and  they  are  valid  in  the  rein* 
forcement  as  well  as  in  the  matrix.  After  adaptation  to  the  boundary  condi¬ 
tions,  the  functions  and  constants  will  have  superscripts  I  and  II  for  fiber 
and  resin  respectively. 

For  the  present  purpose,  equation  (16)  has  to  be  transformed  into 
cylindrical  coordinates  (Appendix  A,  equations  (Al)  through  (A16)). 


BOUNDARY  CONDITIONS  AT  THE  INTERFACE 

A  set  of  conditions  results  from  the  assumption  that,  at  the  fiber- 
matrix  interface,  two  neighboring  particles  of  the  two  materials  remain 


7 


together  during  all  displacements  and  that  the  total  shrinkage  of  the  matrix 
is  greater  than  that  of  the  fiber.  Under  chis  condition,  stresses  will  de¬ 
velop  in  both  materials,  depending  on  the  definition  of  the  displacement. 

It  was  further  assumed  that  the  stress  perpendicular  to  the  interface 
and  the  shear  stresses  at  the  interface  are  continuous. 

There  are  two  possibilities  for  consideration  of  shrinkage.  For 
explanation,  one  cylindrical  fiber  in  a  matrix  cylinder  shall  be  considered 
(Figure  1).  As  a  first  possibility,  it  is  possible  to  describe  the  dis¬ 
placement  with  reference  to  the  free  shrunken  position,  where  one  material 
would  shrink  as  if  the  other  were  not  present;  then  the  displacement  of  a 

point  P  at  the  interface  would  be  in  Figure  1  the  vector  Pjj  P  *  g11 

and  Pj  P  ■  for  material  II  and  I  respectively.  In  this  case,  the 

term  0  in  equation  (9)  should  be  omitted.  The  second  possibility  is  to 

define  the  vector  P0  P  ■  §**  *  51  as  displacement  at  the  Interface. 

The  second  procedure  is  easier  to  use  because  in  the  first  method  the  estab¬ 
lishing  of  directional  signs  for  forces  and  displacements  and  the  two  float¬ 
ing  reference  systems  present  major  problems.  In  equations  (A131)  through 
(A13S),  it  is  shown  that  both  methods  yield  the  same  results. 

In  total,  the  following  boundary  conditions  have  been  assumed  at  the 
interface: 

Displacements: 


sj(a,0,z)  -  sfW.z)  i  -  1,2,3  (17) 


Stresses: 

Oij(a>0>z)  m  Oi.j(fl»0,z)  j  *  1,2,3.  (18) 


GEOMETRIC  BOUNDARY  CONDITIONS  AT  THE  FIBER  HEXAGON 

Referring  to  Figure  2,  it  can  !>e  assumed  for  reasons  of  symmetry  that 
the  hexagonals  will  preserve  their  regular  geometry  during  the  shrinkage 
and  axial  loading  process. 

Therefore,  the  boundary  condition  will  be 


Si1  •  0  »  co80  *  1 


CQ 

COS0 


0  ,  z  I  sin0  ■  e  (19) 


8 


ICOMPOSITE  SHAPE 
•  BEFORE  SHRINKAGE 
P 

°--l 


DEFORMED 


1 


i 


'  I  FREE 


REINFORCEMENT 

CONTRACTION 

FREE  MATRIX 
SHRINKAGE 


(a)  Diagram  Depicting  Components 
Before  and  After  Independent 
Shrinkage 


(b)  Composite  Shape 
Resulting  from 
Combined  Shrinkage 
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(c)  View  A  (Enlarged)  Showing  Combined  Shrinkage  in  Relation 
to  Independent  Shrinkage  of  Both  Materials 


Figure  I.  Reinforcement  Cylinder  Embedded  in 
Hollow  Matrix  Cylinder 
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Figure  2.  Multifiber  Composite  Model 
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where  C  is  half  the  distance  between  two  fibers  before  contraction  and 
o 

e  is  the  total  shrinkage  of  the  composite.  For  comparison  with  the  left- 
hand  side,  the  right-hand  side  in  equation  (19)  has  to  be  represented  by  a 
double  Fourier  series  given  in  detail  in  equations  (A77)  through  (A93) . 

In  order  that  the  fibers  remain  in  the  hexagonal  arrangement,  the  tan¬ 
gential  displacements  at  the  lines  QA  and  OB  have  to  be  zero,  so  that 


§3 ,II(r ,0,z)  -  0  (20) 


^’“(r.J.z)  -  0  .  (21) 


Another  geometric  condition  comes  from  the  fact  that  the  displacements  (and 
stresses)  are  identical  for  every  multiple  of  tt/3  of  the  variable  0  ,  so 
that  in  equation  (5)  the  eigenvalue  must  be  a  multiple  of  6: 


n  ■  0,  6,  12,  18,  etc. 


(22) 


To  satisfy  the  boundary  condition  at  the  hexagon,  which  is  a  function  of  z 
and  0  ,  it  is  necessary  to  have  the  solutions  §  represented  in  orthogonal 
functions  of  z  ind  0  .  Therefore,  we  let  equation  (5)  become 


^k 


k  -  1,  3,  5,  7  .  .  . 


(23) 


The  general  path  of  the  function  5.  is  imposed  by  the  geometry  considera¬ 
tion  of  equations  (A21)  through  (A23) ,  which  express  the  following: 


§1 

is 

even  in  0 

and  even  in  z  ; 

■N 

5a 

is 

odd  in  0 

and  even  in  z  ; 

> 

(24) 

Sa 

is 

even  in  0 

and  odd  in  z  . 

) 

In  order  to  repeat  six  times  in  2tt  ,  the  periodic  characteristic  implies 
that 


^(r.0,2)  ■  §t(r  ,0  +5*  z  j  • 


(25) 
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Since  Che  stresses  are  derived  from  the  displacement  vector  by  derivations 
(equations  9  through  15),  they  also  repeat  six  times  in  2n  period  (equa¬ 
tions  A24  through  A27). 

For  computing  the  shrinkage  alone 
condition  at  the  two  ends  is 


(r,0,±<t) 


The  integration  constants  and  the 
equations  (A44)  through  (A58) . 

The  results  of  the  multifiber  case  show  that  for  the  eigenvalue  n  *  0 
the  stresses  are  constant  with  0  ,  while  for  n  /  0  the  stresses  repeat 
six  times  in  360  degrees. 


SPECIAL  CASE  OF  A  SINGLE-FIBER  MODEL 

In  the  case  of  a  cylindrical  filament  centrally  embedded  in  a  matrix 
cylinder>  the  potentials  expressed  in  equation  (5)  degenerate  into  functions 
independent  of  0  (Reference  7,  page  41-2). 

The  solutions  of  displacements  and  stresses  due  to  differential  con¬ 
tractions  of  the  two  constituents  are  given  in  Appendix  B,  equations  (B8) 
through  (B13)  . 

For  the  axially  loaded  specimen,  the  displacement  and  stresses  are 
given  in  Appendix  B,  equations  (B48)  through  (B53)  .  The  total  stress  and 
the  displacements  in  the  composite  are  obtained  by  superposition  of  values 
from  both  cases. 

The  eigenvalues  for  the  shrinkage  case  are  determined  by  the  boundary 
condition  at  the  two  ends,  where,  for  the  single-fiber  specimen,  equa¬ 
tion  (26)  becomes 


(A40)  without  external  loads,  the  boundary 

“  a“(r,0,±O  »  0  .  (26) 

eigenvalues  for  are  derived  from 


cr33  (c,l) 


ct33  (r  ,1)  -  0 


(27) 


The  only  condition  that  makes  the  solution  for  the  shrinkage  case  nontrivial 
is 


cos  (k-t)  *  0  . 


(28) 
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This  implies  that 


U 

k 


HE 

2 


HE 

21 


n  “  ±1  ,  ±3  ,  ±5  ,  .  .  . 


n  *  1,3,5,... 


(29) 


(30) 


The  eigenvalues  for  the  axially  loaded  single-fiber  composite  are  determined 
from  the  boundary  condition  at  the  free  cylindrical  surface: 


(Til  (b  ,z)  -  0  .  (31) 


The  only  condition  that  makes  the  solution  of  the  loaded  case  nontrivial 
is 


JoM  -  kb  1  ■  0  •  <32) 

The  computed  first  six  eigenvalues  are 

b*ik  -  1.8411  ;  5.3314  ;  8.5363  ;  1.7060 

14.8635  ;  18.0155  ;  21.1643  ;  24.3113  . 


The  boundary  conditions  at  the  interface  of  a  single-fiber  composite  are 
the  same  as  expressed  by  equations  (17)  and  (18)  but  are  not  0  dependent. 


WMERKALRESUITS 

A  single-fiber  model  loaded  at  the  resin  as  shown  in  Figure  3  produces 
axial  stresses  in  the  fiber  which  are  computed  to  be  25  times  higher  than 
the  externally  imposed  stresses. 
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z  « m> 

Figure  3.  Axial  Stress  in  the  Fiber  Produced  by 
Load  at  the  Resin  End  z  =  1.5  inches 
(Single -fiber  composite) 


The  following  values  were  used  in  the  equations  of  Appendix  B: 


E*  =  10  x  10^  psi  (glass  fiber) 

E11  *  0.436  x  10^  psi  (epoxy  matrix) 

v1  *  0.20  (Poisson's  ratio  for  glass) 

v11  *  0.35  (Poisson's  ratio  for  resin) 

a  ■  0.005  inch  (fiber  radius) 

o 

bQ  *  0.6875  inch  (resin  cylinder  radius) 
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The  shear  stress  CTj3  decreases  with  distance  from  the  load  and  be¬ 
comes  zero  at  the  half  length  (z  *  0,  see  Figure  4).  The  axial  shear  stress 
in  the  resin  0*3  (see  Figure  5)  is  a  function  similar  to  a*3  • 

The  radial  stress  in  the  resin  is  shown  in  Figure  6. 

The  tangential  stress  in  the  resin  increases  with  z  and  de¬ 

creases  with  r  (see  Figure  7).  The  tangential  and  radial  stresses  in  the 
fiber  oJa  and  ^n.  respectively  are  equal  and  constant  for  a  certain 
cross  section  and  increase  with  z  (see  Figure  8). 

Shrinkage  stresses  due  to  differential  shrinkage  in  the  same  specimen 
as  above  result  in  apparently  undulated  stress  distribution  (see  Figure  9), 
which  is  more  pronounced  in  a  specimen  with  a  very  high  resin  content.  A 
specimen  with  a  thin  resin  coat  of  a  thickness  of  0.0025  inch  was  also  com¬ 
puted.  Here,  the  computer  results  indicated  no  undulation  of  stresses. 

The  multifiber  problem  equations  given  in  Appendix  A  were  programmed 
for  an  infinite  fiber  bundle.  The  numerical  results  show  that  the  influ¬ 
ence  of  the  solution  with  the  eigenvalue  n  ■  0  ,  Appendix  A,  equations  (A96) 
through  (A99) ,  disappears  at  the  interface  at  a  certain  total  shrinkage  e 
of  the  composite  and  becomes  a  maximum  at  e  “  0  .  It  should  be  pointed  out 
that  e  is  not  identical  with  (3**  ,  the  resin  shrinkage.  If  e  “  0  ,  it 
would  mean  that  the  total  composite  shrinkage  is  prevented  for  some  external 
reasons  (cross  fibers,  fiber  contacts,  etc.).  This  would  produce  very  high 
tension  stresses  in  the  resin.  In  case  the  fibers  touch  each  other,  the 
total  shrinkage  e  would  become  almost  zero  and  the  shrinkage  stress  due  to 
311  would  be  higher. 

The  greater  the  total  shrinkage  e  becomes,  the  more  effective  will 
be  the  solutions  for  n  /  0  .  The  n  /  0  solutions  are  ^-dependent  and 
repeat  themselves  six  times  within  2n  .  The  total  solutions  are,  of  course, 
the  sum  of  both  solutions.  To  indicate  the  influence  of  both  solutions. 
Figure  10  shows  at  the  left  (a)  the  stress  generated  by  the  n  ■  0  solution 
as  a  function  of  the  total  shrinkage  e  • 
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Figure  4.  Axial  Shear  ct13  in  the  Fiber  Pro 
duced  by  Load  at  the  Resin  End 
z  ■  1.5  inches 
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Figure  5.  Axial  Shear  Stress  in  the  Resin 
Produced  by  Axial  Load  at  the 
Resin  End  z  m  1.5  inches 


Figure  6.  Radial  Stress  in  the  Resin  Pro' 
duced  by  Axial  Load  at  the  Resin 
End  z  ■  1.5  inches 
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o  0.4  0.8  l.o  1.2  1.6 

Figure  7.  Tangential  Stress  in  the  Resin 

Produced  by  Load  aJ^CO  at  the 
Resin  End  z  =  1.5  inches 

(single-fiber  composite) 


Figure  8.  Tangential  and  Radial  Stress  in  the  Fiber 
Produced  by  Axial  Load  at  the  Resin  End 
z  -  1.5  inches 

(single-fiber  composite) 
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Figure  10.  Shrinkage  Stresses  in  Radial  Direction 
as  a  Function  of  Total  Shrinkage 


At  the  right  side  (b)  of  Figure  10  is  the  sum  of  the  n  ■  0  and  n  /  0 
stresses  as  a  function  of  0  .  If  the  total  shrinkage  is  not  prevented  and 
the  total  shrinkage  e  assumes  a  certain  value  (in  the  case  of  epoxy  and 
glass,  -0.001132),  then  the  non-zero  solution  is  the  only  one  which  generates 
stresses  at  the  interface.  These  stresses  are  recurring  every  60°.  The 
stress  distributions  around  the  fibers  are  very  sensitive  with  respect  to 
the  total  composite  shrinkage  e  .  An  additional  equation  or  condition  to 
circumvent  or  determine  e  analytically  would  be  very  advantageous.  In  the 
search  for  an  additional  equation,  we  find  that  at  the  hexagonal  boundary, 
the  following  conditions  exist: 


If  the  correct  e  is  chosen,  the  conditions  shown  by  equation  (33)  are 
automatically  satisfied.  The  reason  for  this  is  given  in  the  displacement 
boundary  conditions  of  the  hexagon  (equation  (19)). 

One  can  visualize  that  the  non-zero  solution  is  additionally  influenced 
by  the  distance  2  CQ  of  the  surrounding  fibers  and  the  relation  of  dis¬ 
tance  to  the  fiber  diameter  2  aQ  .  The  closer  the  fibers  come,  the  more 
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pronounced  are  the  circular  six-times-repetitious  stress  distributions 
expressed  by  the  non-ztro  solutions. 

In  Reference  (9),  page  4,  a  condition  similar  to  those  in  equation  (33) 
has  been  used  to  impose  the  hexagonal  geometry  during  deformations.  Using 
this  condition  for  the  present  shrinkage  case,  similar  to  what  was  used  in 
Reference  (9)  for  a  loading  case,  it  would  be  possible  to  set  up  an  in¬ 
homogeneous  boundary  condition  at  the  hexagon  which  would  not  contain  e 
Such  a  boundary  condition  would  contain  only  the  shrinkage  constants  of 
both  materials,  0*  and  0**  . 

The  numerical  results  (see  Figure  11)  show  that  the  radial  stress  at 
the  interface  is  in  compression  at  points  where  the  fibers  are  close  and  in 
tension  where  they  are  farther  apart.  The  shear  stress  (Figure  12)  at  the 
interface  is  a  maximum  at  $  “  22s  and  zero  at  0  ■  O’  and  30°  .  The  total 
shrinkage  displacement  between  two  fiber  axes  in  the  model  was  2  x  (-0.001132). 
The  dimensions  of  this  particular  multifiber  specimen  were  aQ  ■  0.14  inch 
and  C0  ■  0.20  inch  . 

Appendix  C  describes  the  methods  used  and  the  results  of  a  photoelastic 
determination  of  stresses  of  selected  single-fiber  and  the  multifiber  models. 
They  are  represented  as  stress  differences.  It  is  now  possible  to  present 
the  theoretical  analysis  in  the  form  of  principal  stresses,  rather  than  in 
the  form  of  stress  differences,  and  to  compare  the  theory  with  the  experi¬ 
mental  values. 
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Figure  12.  Shear  Stress  at  Interface  Due  to 
Shrinkage  aQ  =  0.14  inch;  CD  = 
0.2  inch  (Components:  glass  and 
epo 
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APPENDIX  A 


TRIDIMENSIONAL  ANALYSIS  OF  THE 
UNIDIRECTIONAL  MULTIFIBER  COMPOSITE 


LONG  CYLINDER  WITHOUT  AXISYMMETRY 

The  displacement  and  stress  components,  in  terms  of  Papkovitch  func¬ 
tions,  in  rectangular  coordinates,  are  quoted  as  follows: 


_L_2_ 

4-U-P)  IK 


(XaYas*  ?.) 


(AI) 


To  convert  these  expressions  into  a  system  of  cylindrical  coordinates, 
one  can  perform  the  following  proper  and  admissible  transformation. 

One  transforms  from  (x  ,  y  ,  z)  to  (r  ,  0  ,  z) 

e°  (r>,  r*,& 

(u,  v,w)  to  %) 

%%%  1  <&,  w 

The  results  are  collected  as  follows: 
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for  displacements, 


f ,  “  P,  “  +if!  +  ?.)J 


for  stresses, 


a 


'a 


4.  <3-.  r  0v7*  ^  i  ci  _  _j _ 

r m-p)L^  vc  ^c^hk  fisi^ 

•(*?,+*£+  fc) 


+2(M)[^v*e  “(5^^:  +  a^^/] 

•f/^F4^  +fJ  (A7) 
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for  Laplace's  equation  of  Papkovitch  functions, 


v;  (f(  <^4?  — — 0 

+J>,  =  6 

Vj  fi-0 

V‘  ?.  =  o- 


(A12) 

(A13) 

(A1A) 

(A15) 


A  solution  of  Laplace's  equation  in  the  finite  cylinder  in  the  domain 
(0  <  r  £  a  ,  0  £  $  £  2tt  ,  0  £  z  i  ()  takes  the  following  general  form: 


■f"  (^#0  + 


(A16) 


i 
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Combination  of  equations  (A12)  to  (A16)  gives  the  folloving  results  of 
Papkovitch  functions: 


-f F„,  csafaVj  <u s4 

'tmi 

mSin(*>tO  +  Coj (w^) J 

■f  C  C..  i  D*,  *)  (  An  +  kj-*) 

t  (  l"t>  f  fiL  *)  (G*>f  hl»  tjk)  St*j> 

+  Tlf  [A*  JWUA)  1 6*  r^4] 

*[C*  s4  OWm,*)  4 

E«( Sin  (oi+i)HSih(to-/)4>) 

iFJ  ('n+O'H  C*J('H)  49] 
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(A17) 


~~  1  jZ  [a^^a) + #*,  x^a)J 

‘  [  C„Aj  Si  m  +J^j  Ce}U^^)j 

*■£[ 

4Fn(s,v^)4-siVW^)J 

+  T  ( f  [  (^VA*  4)  A  ^  YJ^  jo] 

SYh*c  l  « 

•  [  Aim  (Kn  ^ J 

'  i[A/w(SiV'»>$<Ms»Vi  ('’”-/)$) 
+  Pw  ( Ccjf'Wiwty-f  eofw'O^J 
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*[a/w  Si‘*i('»l4>)  +  J5wC<5f^+)J 


'£  B„  S1V1  (01^)  +  p„  tcj  (<w<£)J  siw^ 


f  (  Let  ~t~  Meo  ^)((ree  ■f  Wm  /»4-k)  C«4 
—  (Cec,  *f  JPoe  *)  ("  A»e  +  6ee  A)  SlM  ^ 

(J  (A18) 

d 

ft  =T{l  [(V4V+M(*^J 

.[VW^  ^  c®Us^,*jj 
*  [  S?"  (5  Xj^ 

e^o 

-f  2  [  &,Aa  f  ^~l]  (  Ss°+  Vs] 

* [  M^Sim  (5.4>)  f  Xj_ ^cd^Oj 

f(5*e+ T^X <&*+&,  ^A.)  (A19) 
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Domain  and  Boundary  Conditi.ns 


Domain: 


(0  <  A-  ^  00 


< 


Ck  ^  A.  ^  <Vc«3^ 

o  ^  4  <  % 


^  < X 


for  fiber 
for  resin 


(A20) 


Geometry  consideration: 


For 


( A.^, —  S;( ( £.,-<}> 2) 

(  A-/  ^  ^  (Ay  4,  <A21> 


[  t^K  *,  *)*=*“£.( A-ri  a) 

1  a,  4»  *)  =  §*  -*) 

(A22) 

f  i3  ( &/ 4>,  ~ ^ 

§3  (Ay  4>  S>)  ==■  §3( Ay- 1 *v> 

J  (A23) 
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Periodic  characteristic  of  tt/3  period  implies  that 


4>,  *)  =  f  l  [fL  ,4>  +  *),  CA24) 

^  (*,  4»  a)  =  ||  ( . a,  ++ £  *> ;  <*«> 

and  QTx  (  ^  4>  -2>)  =  *>),  (A26) 

A,  =  fpN  M> + *  *)•  <«» 

Displacements 

Su^t  Ltution  of  equations  (A21)  to  (A23)  into  equations  (A3)  to 
(A6),  (A9) ,  and  (All)  and  taking  into  account  the  consideration  of 
the  previous  section  yields: 


Radial  displacement  for  fiber, 


33 


34 


L 


+r  ■>*,]  ft-®’--"*"1 

■*+  ^  ^ 

**  a/h + (/a  '7tA'"  j 

'  HO 


-(Cl  *‘a‘'  ■KAiia-^wi*.-': 


(A29) 


Tangential  displacement  for  fiber, 


•A”+l  z*]5.V^) ; 

v  me?  -J  mm 
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ii — 1^4^ 


Tangential  displacement  for  resin, 


l(A*)we  ( *+4(i-jf))  A™ 

6Uh,  '71  ^ zl +  (<£*)„,'*  A~* 


(A31) 
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+r[(o>- 


4P*)h! 


Cjryty^>) 


(A32) 
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Cx=?  t  oo 


Xl  |f 

+rif(  cl%^  xm^« 

"7,sf^mUHitn  ‘  /  /  ^ 

>?  /■  60  rv 

f[  [  (Q 

Mr  lty$«  '  '  ^ 

^IJT  (C]/3  ^ 


C*Q  r  ^o 


^4* 


(-/u 


?)r*C<56?xh) 


■+Z  f[  (KL(y^ 


i  [  f  (<L  f*~* + WIltM**)  A-”] 

-w  J 


(  %)'  CcS^^fc) 


(A33) 
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Stresses 


Radial  stress  for  fiber, 


40 


.XSimU'W**) 


JkO 

+  >  (A*)  Wl-^-^-D^coi (<*4) 


'D*i>ly" 


Radial  stress  I  or  resin  has  all  terms  defined  in  the  same  form  as  equa- 
t '  on  (A34)  except  for  replacing  all  superscripts  I  by  II,  plus  the 
t  •  rms  involving.  Bessel  functions  of  the  second  kind  and 

J  (  ft1]  fa 

£—  *lyHO  Jtm 


Shear  stresses, 
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rr  i 

(A* JM0  "» (' a*'1  l s'iM ■ 


(A35) 

_ 

has  all  terms  defined  in  the  same  form  as  equation  (A35)  except  for 
replacing  all  superscripts  I  by  II,  plus  the  terms  involving -Bessel 
functions  of  the  second  kind  and 

h =£/" 
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f  P-o 


•si 'Cnfafy 

+J  ■[  (  Cl 

-  -77  f 

~^-nAj'Jh-ij/^uijft-) 

+  '&C*jp0 

■h 

003(9?  4>) 


(A3  6) 


G^has  all  terms  defined  in  the  same  form  as  equation  (A36)  except  for 

replacing  all  superscripts  I  by  II,  plus  the  terms  involving  Bessel 
functions  of  the  second  kind. 

Axial  stress  for  fiber, 

<sT-  £[f {f 

#  CX55 1 1  \  •  C<3  (9i 


r  bo 


+Z]Z 
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•Coitanty) 


6-0 

+  ^  (Al\o  h"  c<5 ('nty 

/H/V" 

■7)=U>/i/" 


(A37) 


Axial  stress  for  resin, 


has  all 
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terms  defined  in  the  same  manner  as  equation  (A38)  except  for  replacing 
all  superscripts  I  by  II,  plus  the  terms  involving  Bessel  functions  of 
one  second  kind  and 


Boundary  Conditions  at  the  Interface 
The  continuity  of  displacements  implies 


$,*<«,  4,*) 

The  continuity  of  stresses  implies 


(A3  8) 


(A39) 
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Axial  stresses  at  the  ends  should  vanish;  i.e., 

^(K  <*/*)  -  <£(jlA,  W~°- 


The  sides  of  the  hexagon  should  remain  fixed;  i.e.,  . 


&/  4/  ^  “ 5*  +/ $si*+  ~e  ■ 

(A4i; 

Since  is  approximately  57.  of  ,  we  then  get 

£i*(cis$,  t £)  =  6  <M2) 


or 


(A43) 
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Equations  to  Determine  Coefficients 


The  constants  to  be  determined  are  listed  as  follows: 


f-»  =  O,  L,  /2y  "  ' 
{  A.-/,  3/  Sy 


(  -»  4  !Z/  /g>  "  ' 


) 


These  constants  will  be  determined  by  the  following  equations. 
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Combination  of  equations  (A30),  (A31),  and  (A38)  yields 


-  ( AI) J>  ++f,~^&  ^ 

—  Y hh^Ah**/1)  - 

Combination  of  equations  (A32) ,  (A33) ,  and  (A38)  gives 

(A  1/)^  (rfi 

+  (AU*(yu,^^<l) 

t  ( A 

+  ($)*  ( /V*a)  Yn-,y*oH*x) 

t  (  %L(-/l»4  Y*,yu*)  =  *  • 


0  .  (A44) 


(A45) 
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/~jp* 


I 


I 


■+/^«h)-4  ^  Yn-ij^th-oi/ 


t-ZP*  ^ 


■+ ^  Y*+i 

.  j^W) -(vHX'M'Khj?)  Xs<*$ 
l-zf  J 

(A46) 


Equations  (A35)  and  (A39)  yield 


( A*L(iS)  ^ ^4-^A) 

7*«W  * 

'Jff 
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i 


51 


52 
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Equations  (A34)  and  (A39)  give 

t  <r"'*( ;£)  -o 

(A52) 


(a;i 


-2 

06  y^x 


+ (€), 


,  E^r 

00  l+P* 


~  (AojJoo  A  ji^  ^  (^jJoo  ^  /^*  ^/i-  HP* 


_  E 


l-zjr  f  /-5^P  • 

Equations  (A35)  and  (A39)  yield 


M‘i.  (7^5)  ■»  ( '»-tu-P‘»  aM 

+  («£)„(& 


(A53) 


(A54) 
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Equations  (A37)  and  (A39)  give 


(K)cc  4^(1- IP*)  +  Pi  z(  '-^) = 

and  (A55) 

(A^0  = -o+p^f- 

Equations  (A30) ,  (A31) ,  (A39) ,  and  (A40)  give  A56) 


~KL  ^  ~*{,~pJL))  *'H  ~ 


(A57) 


Equations  (A29)  and  (A42)  yield 


55 


(A59) 


«L*o-»r)(£f)+(eZX.(}£) 
=  . 
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Transformation  of  Coordinate  Systems 


One  has  the  following  transformation: 


o  3,  p — - ( }) 


57 


A.  ~JX  +  f 

4®^  '($) 


(A61) 


(A62) 


(A63) 


(A  64) 
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From  equations  (A03)  to  (A66) ,  we  get 

■&x  24.  +  VC1&  t  Vc  ^ 

—  C<ts<i>^ _ -Sifljfe  2L 

*ta  a  ?>4> 


3J- 


■-  S'«*\  + 


1-  =1^:  +  t>4>  g,  +^1_ 

2%  ^%2/t  +  +  >*<)* 

=  ^_. 

7>* 


(A68) 


(A69) 


The  relationships  between  two  coordinate  systems  for  displacements  and 
Papkovitch  functions  are 


U  =  §,  e.«4>  ~  S23in$ 

V  —  §t  <sinty  + 

^=13 


(A70) 


(A71) 
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Combination  of  equations  (A3),  (A70) ,  and  (A71),  and  dropping  the  P 
term  yields 


U  —  Si  C0S4- 

~  f  cas<t>u  [^4* ( T'^- 

iKS'i*$(fi  s",4>  +  p» 

S-k4>) 

* 

+4s.k4>  f|>l5.'4+f^‘l>)+*^]J 


CC3  <\>  —  ^  5fW 


Similarly,  combination  of  equations  (A4) ,  (A70) ,  and  (A71),  and  dropping 
the  PQ  term  yields 
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V  “  I,  c<o4- 


pi  Si*k4>  + 

-OhsI^+Mw+zk] 

+c^fAp+*r»]j. 


(A73) 


nui  dp  lying  (A72)  by  co«0  and  (A73)  by  .In*  ,  adding  togel 
thinking  use  of  the  trti«u»etric  Identic/  sin®* 


5,  ~  Pi  fikty  f fe  fA-f,  t  •  «»*) 

Similarly, 

4.  ““  J*  ~  J  A-P* +  *J?j[  <*”> 

^  js  “  J(0~P)  J^.[A-fi  (AW) 
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Fourier  Representation  of  Boundary  Condition  Equation  at  the  Outer 
Surface  of  the  Hexagon 


At  the  outer  surface  of  the  hexagon,  one  has 

or 

5  (  ^Ty+/*)  =¥ 

We  want  to  represent 


(A77) 


(A78) 


5  .i-'-c  $ 


''Yh^O 


+  Z 


~7V*-( 

6o 


+  Z 


<'7> f=») 

+rr^ 

•'hh  Ah 


Now  define  a  function 


(A79) 


(A80) 
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The  coefficients  of  the  Fourier  representation  are  as  follows: 

0  ^ 


o 


^  fjtsin  [(i-m-t)  %]  _  X  Sin  [(C*~3)%] 
l  Cin-i  £*»-j 

-f  _  f*  cgjWM  .fik 

(>y*S  >Ja 

“  |j±T  -  fife  +  <±? 

=  £  fffirfess  ■ 


(A8S) 


When  Ms  j, 


(A86) 


When  m  Z  , 


If*  *?SP*“  Mr  -/  *s2Hf 


“^f(//«?7  '  <f-3l ' 


(A87) 


Continuing  to  attune  nunbert,  one  gett 
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‘#U© 


=  «r 

7T  l 


s 


8 


^  llm  -tfUm-ty 


+  ( 


|^(|  Tu-lai-mj  +  H*i  V] 


»h«n 


(AM) 


whan  -1H  «  Ctfe*/ 


v3)<«, 

rr 


^  </*J  «t>  ii^  J*J+ 


(A89) 

(A90) 
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when  £  =  odd 

when  even 

(A91) 


4)o(mAJ 

yj"  &  Ate  ^  <A<^> 


when  “7^7  =:  ei/e*.  . 
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PLANE  STRESS  FIELDS  WITHOUT  POLAR  SYMMETRY 


The  solution  to  the  biharmonic  function  in  cylindrical  coordinates, 
for  the  stress  fields  without  polar  symmetry,  is 


L—£  (A^i 

/t^'"  (A94) 

^  As  ton  h.  +AiA*+A7$- 
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/ 


where  F(r)'e 
section. 

ere  functions  of  r  only  end  ere  defined 

in  the  next 

Definitions  of 

—( 'h-b'h  A*"4 

(A100) 

& 

II 

C-w-zt)  A* 

(A101) 

F„»)  - 

--7J  (/?»+()  A.-*''1 

(A102) 

F.W- 

h.~K 

(A103) 

F*ui— 

( <*>  -O'n  a"’* 

CA104) 

(-71  ■+&('»+*)  /h~ 

(A105) 

£ 

II 

■  m(n HOti*'* 

(A106) 

Fsn  (A)  — 

■■  (Oi'Z)fo-i)  A’*' 

(A107) 
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FJJi)  *=('n-0-v  kn~*' 

(A108) 

FJ*  - 

'»('*+!)  JlK 

(A109) 

F-W  = 

-Mortd 

(A110) 

— ( “W  -/)-»)  yt""' 

(Alii) 

II 

_  lllHjlDy-i 

(AU2) 

1 

£ 

1 

i 

?*• 

* 

(A113) 

^,(<0== 

22ii±i^  a -«  -/ 

E 

(A114) 

F/nfft)  = 

£ 

(A115) 

F/7, '*•)=- 

*2i£l±£lh»-i 

E  11 

(A116) 

'»(l+J>)+*t  v«+» 

E 

(A117) 

E  A 

(A118) 

5«  fo)  = 

-7)  (/■+#)-#  l-H+I 

E 

(A119) 
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1.  At  the  interface,  the  stresses  and  displacements  should 
be  continuous;  i.e., 


A«  F>) 

-a '.FZM-AT.F’U-O 


x  _  x 


'jn'  3/t 


(A120) 


A*F*M  fAl  F,» -AlKjO-AFKj*) 


x  — x 


'/*  *  <?« 


-A* 

/x3*'  Hn  4*  1  H*  ' 

•  I1tk4>)=^(a,4>; 


J!*  '  IOH. 

A  X-  r-X 


in  qn 


JLt t' 


a  r.  f;^ +a^-a:o-)-ax^> 


(A121) 


(A122) 
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40  =£*K4) 

-  4*  F<t (0i)  ~K  £ ^  ^ 1 5  •  (/ 


2.  At  the  outer  surface  of  resin,  it  should  remain  a  regular  hexagon 
after  deformation;  i.e., 


( ^,4)^4 


FZ'(£f)+AlF*(£f) 
rff  -(-0^8  ,f  fj, , 

when 


tt[|  (u-d;*)(ii4H'd 


(A124) 
when  ^Th^JUjJAt 


'M^f,  2,  3,  •" 
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Boundary  Conditions  for  n  ■  0 

1.  At  the  Interface,  the  stresses  and  displacements  should 
continue . 


*•  o)  =  <y*(OyO)  . 


It  implies 

4X= 0 

(A126) 

2  A.  ~AfX <*~z  —  • 

(A127) 

b-  *  <$*(4,4)  ■ 

It  tallies 

II 

II 

«  - 

CA128) 

«•  $*(^4*)  ~  ■=0  • 

This  yields 

+  Af*^^:'  =  6. 

(A129) 

2.  At  the  outer  surface,  the  hexagon  remains  regular  after 
deformation. 
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3.  Comparison  of  Two  Reference  Systems 


With  reference  to  the  original  position  of  a  particle  before  shrink¬ 
age,  the  shrinkage  stress  in  the  shrunken,  equilibrium  position 
is,  from  equation  (9), 


(A131) 


r  *  a  =  composite  interface 
before  shrinkage 


Figure  13.  Graphs  Showing  the  Comparison  of 
TVo  Reference  Systems 
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With  reference  to  the  location  where  the  material  particles  would 
shrink  if  the  other  material  would  not  be  present,  the  tensor  equa¬ 
tion  (9)  for  stresses  in  the  matrix  in  one  dimension  becomes 


(A132) 


In  Figure  13,  the  relationship  between  the  displacements  is 

-X* 

r  -  f-  s 

'dn  ' 


Substituting  equation  (A3)  into  equation  (A2), 


(A134) 


shows  that  both  reference  systems  yield  the  same  shrinkage  stresses 
an  ,  equations  (Al)  and  (All).  The  comparison  could  be  made  in  three 
dimensions  and  would  show  that,  in  general. 


■  <ai35) 
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APPENDIX  B 


IMPIMEHSIWAL  ANALYSIS  9F  IH£ 
SMgLE-fIBER  CfflgflSHE 


SH&MKME  CASE 

The  solutions  for  a  biharmonic  equation  in  a  system  of  cylindrical 
coordinates  with  axisymmetry  are  as  follows: 

L=z  i 


where  I0(4kr)  and  Ix (n^r)  are  modified  Bessel  functions  of  the  first 
kind  and  of  the  zeroth  and  first  order  respectively;  and  and 

are  those  of  the  second  kind  in  their  corresponding  zeroth 
and  first  order. 


Pis  placement  8.  and  Stmm 


The  displacements  and  stresses  in  terms  of  strain 

function  are: 

E  =  -!±2>XL 

(B2) 

°i  E 

(B3) 

(BA) 

<£-&(»*!■ -tik) 

(B5) 
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(B6) 


Combination  of  equations  (Bl)  to  (B7)  yields 


(B13) 


where  (r)  ,  .  .  .  ,  Fis  (r)  are  functions  of  r  involving  IQ(p,  r)  and 
Ii(ujcr)>  and  Gi(r)  ,  .  .  .  ,  G1^(r)  are  functions  involving  K  r)  and 
Mukr).  ° 
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(B14) 


Definitions  of  F(r)  and  G(r)  : 

Fia/V  =  ~  ^ 

F#fA)  =  ~(^9  ^  <B15) 

FA  =  (^A  (B16) 

FJa)  =  [4  j-  Uit  m]  (bw) 

-  A  (BIS) 

f>/A)  =/A  [  f -^)  J, ^a)  <m9) 

FjK^A)  ~  ^/K(J  (B20) 

faff)  ~  f'u  (2J)~i)  Xjfaty  (B2 1  > 

Ffl/d  —  A^  (B22) 

FA  =/A[*  ^-P)  X -AH  Aj[  (B23) 

RifcfO  -/^  JiMA)  <B24) 

F*A  [A  4  t2(M;j^wi[  <b25> 

G,Uiifi=  K|(A^)  (B26) 
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Superscripts  I,  II  were  omitted  in  equations  (Bl)  to  (B37). 
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Imposing  Boundary  Conditions 

1.  At  the  interface,  the  continuities  of  stresses  and  displacements 
are  assured. 


k. 


b)  I31('«/*)  +  §/('«;,*)  =  %>(/$*- (3*) 


0  q;xK*)  -  <5^*) 
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4)  ®/3  (*/*$  ^3 

A*  FM' +  AM> F*,®  '  -fifCS*) 

2.  At  the  outer  surface  of  the  matrix,  the  resin  should  be  free  from 
normal  and  shear  stresses. 

a)  o;xa >,z)=*o 


(B42) 

t) 


(B43) 


3.  At  the  ends,  axial  stresses  for  both  fiber  and  resin  should 
vanish . 


0^7  A,  jO  =  =0 


The  only  condition  that  makes  the  solution  nontrivial  is 
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(B44) 


*=  ±t,  ±-5,i£"> 

(B45) 

a.—  f  /  3/  v  v  »• 

(B46) 


Remarks 

1.  In  deriving  equations  (B38)  and  (B39) ,  we  have  used  Fourier 
expansion  for  a  constant  and  a  function  z  . 

2.  It  is  interesting  to  notice  that  the  eigenvalues  determined 
by  equation  (B46)  are  consistent  with  the  argument  of  the 
orthonormal  function  used  in  the  Fourier  representations  in 
equations  (B38)  and  (B39) .  The  coincidence  is  not  accidental; 
rather,  it  is  intentional. 

3.  The  constants  to  be  determined  for  each  k  are  A*  ,  A?  , 

II  TT  II  II  *  J 

Ai  ,  A£A  ,  A*  »  and  A£  .  These  six  constants  can  be  found 

by  solving  the  simultaneous  equations  (B38)  to  (B43). 
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AXIAL  LOADING  CASE 


The  superposition-product  solutions  of  the  biharmonic  equation  in 
cylindrical  coordinates  with  axisynxnetry  has  the  following  general  form: 


(B47) 


where  an<*  ^i(^kr)  are  cylindrical  harmonics  of  the  first 

kind  and  of  the  zeroth  and  first  order  respectively,  and  and 

Yi(^kr)  are  that  of  the  second  kind. 


Stresses  and  Displacements 
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"I 


(B53) 


where  Fx  (r)  ....  .  Fia  (r)  are  functions  of  r  involving  J0(jxkr) 
and  •Jj/nic1)*  and  Gj^  (r)  ,  ,  Gia  (r)  are  functions  of  r  involving 

YoKr)  and  YiKr)  • 

Definitions  of  F(r)  and  G(r)  : 
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(B61) 


F^ja.)  ==  Xi>^ /k3 

=  -/iC'J.fr*')  <»«> 


FiJfi  ~/£  {■*(*-■»)'%  (/<aA)  yi^y<tA)l 

FJ*)= /tiffin*) 


(B63) 

(B64) 


(B65) 


Gx  (r)  ,  .  .  .  ,  (r)  are  defined  exactly  in  the  same  manner  as 

equations  (B54)  to  (B65)  except  replacing  J^nkr)  and  Jx  (p^r)  by 
Ko  (m-R17-)  and  Kj/p^rJ  respectively.  Superscripts  I  and  II  were  omitted 
in  equations  (B47)  to  (B65) . 


Impgsins  Bwtrtary  tanlUlon? 

1.  At  the  interface,  the  continuities  of  displacements  exist. 
This  implies 
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b)  §/(<*,*)  =“§/K*) 

Aki  fwa) + ^ Fn/*) ~a  Ft  ^  ~\kF^  ~ 0 

2.  At  the  outer  surface,  radial  stress  of  resin  should  vanish; 
i  •  e  • , 

o;,  fty  as)  =o  . 

It  implies 

At'Fj/W  +  • 

3.  At  the  ends,  resin  is  subjected  to  axial  loading;  i.e., 

c&Va-,  0  =  Fa;  —  load  per  unit  area  . 

This  yields 

■/a  ^ [A i *  - A^u -£>] - q, 

and 


(B67) 


(B68) 


(B69) 


(B70) 
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(B  71) 


Equation  (B70)  Implies 


Then  combination  of  equation  (B69)  and  (B71)  gives 


(»72) 


Determination  of  Eigenvalues  and  the  Dial -Bessel  Representation  of 
Loading  at  Ends 

1.  Eigenvalues  are  determined  from  equation  (B68) ;  i.e., 


2.  The  loading  at  ends  can  be  represented  as  Dini-Bessel 
expansion;  namely. 


(B73) 


(B74) 


=2 

A 

where 


(B75> 
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(B77) 


Remarks 

1.  Eigenvalues  ^  (>  0)  are  determined  by  equation  (B74) . 

2.  Coefficients  of  Dini-Bessel  expansion  for  loading  at  ends 
can  be  found  from  equation  (B76)  or  (B77). 

3.  Constants  to  be  determined  are  ,  and  A*£  .  These 

values  can  be  found  by  solving  equations  (B66) ,  (B67),  and 
(B72)  simultaneously. 

4.  The  condition  that  a/b  is  much  smaller  than  unity  is  assuned. 


APPENDIX  C 


PHOT  OE LAST  IC  EXPERIMENTS 

SINGLE-FIBER  SPECIMEN  PREPARATION 

To  perform  an  experimental  evaluation  of  the  shrinkage  stresses 
for  a  single  fiber  and  epoxy  resin  matrix  via  photoelasticity,  four 
cylindrical  specimens  were  fabricated.  Typical  dimensions  of  the  cylin¬ 
ders  are  shown  in  Figure  14.  Specimen  one  contained  no  fiber;  specimen 
two  contained  a  4-mil-diameter  boron  fiber;  specimen  three  contained  a 
5-mil  glass  fiber;  and  specimen  four  contained  a  10-mil  glass  fiber. 


With  Dye thynol amine  (15  phr) 

Used  as  Resin  Matrix 

Figure  14.  Typical  Dimensions  of  Cylindrical  Specimens 

The  method  of  fabrication  for  the  cylinders  and  fibers  was  to  locate 
a  fiber  in  the  center  of  a  metal  tube  (the  fiber  was  supported  to  allow 
axial  freedom  of  movement  as  the  resin  cured)  and  to  pour  Shell  Epon  815 
resin  with  dyethynolamine  15  PHR  (parts  per  hundred  weight  ratio)  slowly 
in  to  fill  up  the  tube,  allowing  air  bubbles  to  be  forced  to  the  sur¬ 
face.  A  rubber  stopper  was  used  to  plug  the  bottom  tube  end  (the  fiber 
end  simply  "rested"  on  this  stopper). 

The  resin  cure  temperature  cycle  was  200°F  for  12  hours,  with  post¬ 
cure  at  250°F  for  4  hours.  A  parting  agent  applied  to  the  inner  surface 
of  the  tube  allowed  the  resin  to  shrink  during  the  curing  process  and 
during  the  cooling  period. 

Figures  15  through  18  show  the  specimens  photographed  under  direct 
polarised  light  (quarter  wave  plates  are  in  the  system  to  allow  only 
stress  magnitudes  to  be  seen).  White  light  was  used  for  the  photoelastic 
analysis  because  the  various  colors  indicate  very  clearly  the  stress  mag¬ 
nitudes  present  in  the  specimens . 
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Figure  15.  Epoxy  Resin  Cylinder  without  Fiber 
Viewed  under  Polarized  Light 
(Quarter  wave  plates  are  in  the 
optic  system) 
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Figure  16.  Boron  Fiber  (0.004-inch  Diameter) 
in  an  Epoxy  Resin  Cylinder  Viewed 
under  Polarized  Light  (Quarter 
wave  plates  are  in  the  optic 
system) 
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Figure  17.  Glass  Fiber  (0.005-inch  Diameter) 
in  an  Epoxy  Resin  Cylinder  Viewed 
under  Polarized  Light  (Quarter 
wave  plates  are  in  the  optic 
system) 


93 


Figure  18.  Glass  Fiber  (0.010-inch  Diameter) 
in  an  Epoxy  Resin  Cylinder  Viewed 
under  Polarized  Light  (Quarter 
wave  plates  are  in  the  optic 
system) 
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As  can  be  seen  from  Figures  15  through  18,  the  boron  and  5-mil 
glass  fiber  buckled  due  to  the  compressive  load  applied  to  the  fiber 
by  resin  shrinkage.  As  nearly  as  can  be  seen,  the  10-mil  glass  fiber 
specimen  did  not  buckle,  but  remained  straight. 


From  the  photographs  of  the  boron  fiber  and  5-mil  glass  fiber 
specimens,  the  following  information  was  obtained: 


Boron  Fiber  (4-mil  diameter) 

Buckle  wavelength  ■  0.238  inch/cycle 
Double  amplitude  ■  0.0125  inch 

5-mil  Glass  Fiber 

Buckle  wavelength  ■  0.212  inch/cycle 
Double  amplitude  B  0.015  inch 


Since  the  theoretical  analysis  is  predicated  upon  the  fiber's  not  buck¬ 
ling,  the  specimen  chosen  for  quantitative  stress  evaluation  was  the 
10-mi-l  glass  fiber  specimen. 


In  order  to  photoelastically  analyze  the  stresses  in  the  resin 
around  the  fiber,  a  thin  slice  along  the  fiber  axis  from  a  cylindrical 
specimen  must  be  used.  This  is  to  allow  the  polarized  light  to  incident 
perpendicular  to  the  specimen  surface  and  also  to  make  the  stress  field 
appear  two-dimensional,  by  having  the  thickness  of  the  specimen  small 
enough  that  the  thickness  stress  component  can  be  neglected.  To  cut 
slices  for  analysis  from  a  specimen  at  points  and  in  directions  of  de¬ 
sired  stresses  is  a  standard  technique  for  photoelastic  measurements. 

A  problem  is  encountered  with  slicing  when  thermal  stresses  are  the 
quantities  to  be  measured,  since  during  cutting,  heat  is  produced  and 
this  perhaps  distorts  the  original  stress  values.  In  addition,  the 
remainder  of  the  resin  material  might  result  in  stress  releasing  of  the 
fiber. 


It  was  therefore  decided  to  slice  the  10-mil  fiber  cylinder  in 
three  steps  in  order  to  observe  this  effect.  The  cutting  was  also  done 
under  a  high-pressure  water  jet  to  minimize  the  effects  on  the  original 
thermal  stresses  present  in  the  specimen. 

It  was  decided  to  cut  the  specimen  first  to  1  inch  across  the  flats, 
then  to  1/2  inch  thickness,  and  finally  to  1/8  inch  as  a  goal  for  analy¬ 
sis  evaluation.  A  diagram  of  the  cuts  is  shown  in  Figure  19. 
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Figure  19.  Typical  Cuts  Made  on  10-mil  Glass  Fiber  Specimen 
to  Obtain  a  Specimen  for  Quantitative  Photo¬ 
elastic  Stress  Evaluation 


The  stress  patterns  of  these  cuts  as  shown  under  polarized  light 
(quarter  wave  plates  in  the  system)  are  found  in  Figures  20  through  22. 
Additional  slices  cut  parallel  to  the  fiber  but  perpendicular  to  the 
first  cut  mi^nt  have  further  reduced  the  effect  of  heat  during  the  cut¬ 
ting  operation,  however,  this  idea  was  rejected  because  such  slices 
would  have  cut  off  the  zero  Isochromatics  reference  points  which  are 
necessary  for  the  quantitative  analysis. 

The  pictures,  of  course,  show  some  change  in  stress  pattern  between 
the  first  cut  and  the  final  cut.  One  has  also  to  be  very  careful  in 
evaluation  of  the  fringe  lines  and  remember  that  the  fringe  order  in¬ 
creases  with  thickness  of  the  specimen. 

An  elementary  examination  revealed  that  the  stress  near  the  center 
of  the  fiber  for  the  first  cut  and  the  stress  for  the  final  cut  turned 
out  to  be  fairly  close  to  each  other  in  magnitude.  On  this  basis,  it 
was  concluded  that  the  original  stress  pattern  was  not  greatly  harmed 
by  slicing. 

It  was  decided  to  use  the  Budd  Company  reflection  polarscope  as 
shown  in  operation  in  Figure  23  to  evaluate  quantitatively  the  stresses 
in  the  slice  from  the  cylindrical  10-mil  glass  fiber  specimen. 
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Figure  20.  Glass  Fiber  (0.010-inch  Diameter) 
in  Epoxy  Viewed  under  Polarized 
Light  after  Slicing  the  Cylinder 
to  1  inch  Across  Flat  Areas 
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Figure  21.  Glass  Fiber  (0.010-i:  ch  Diameter) 
in  Epoxy  Viewed  under  Polarized 
Light  after  Slicing  the  Cylinder 
to  1/2  inch  Across  Flat  Areas 
(Second  cut) 
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Figure  22.  Glass  Fiber  (0.010-inch  Diameter) 
in  Epoxy  Viewed  under  Polarized 
Light  (Direct  polarscope)  where 
Distance  Across  Flats  is  1/8  inch 
(Third  cut) 
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Figure  23.  Budd  Company  Reflection 
Polarscope 
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EXPERIMENTAL  ANALYSIS 


The  reflection  polarscope  operates  by  observing  the  fringes  in  a 
birefringence  model  coated  on  one  side  with  a  reflective  coating.  In 
effect,  the  reflection  polarscope  simply  doubles  the  specimen  thickness 
since  the  light-path  thickness  is  doubled. 

The  first  step  In  quantitative  stress  evaluation  by  photoelasticity 
is  to  determine  the  stress-optic  coefficiency  for  the  particular  material 
of  the  model. 

This  is  done  by  using  a  standard  calculable  stress  distribution 
model  and  observing  the  fringes  as  this  model  is  being  loaded  (this 
technique  is  discussed  in  detail  in  Reference  1). 

For  this  case,  a  flat  sheet  of  the  same  thickness  as  the  final  cut 
on  the  10-mil  glass  fiber  specimen  was  cast  from  the  same  resin  and  with 
the  same  cure  cycle.  A  "dogbone"  tensile  specimen  was  then  machined 
(under  cooling)  from  this  plate. 

The  test  section  of  the  machined  dogbone  was  coated  with  an  aluminum 
paint  for  reflection  of  the  white  light  used  in  the  Budd  Company  instru¬ 
ment. 


The  dogbone  specimen  was  then  placed  in  a  tension  testing  machine 
and  load  was  applied.  As  the  load  was  Increased,  the  "tint  of  passage" 
or  fringe  order  was  recorded.  Several  load  and  unload  cycles  were  used 
in  this  process.  In  every  case,  the  applied  stress  was  below  the  elastic 
limit  and  conformed  to  Hooke's  law. 

A  plot  of  the  fringe  order  versus  applied  load  gave  the  slope  of 
the  load-fringe  order  curve.  The  stress-optic  coefficient  was  then  com¬ 
puted  as  shown  in  Reference  1. 

In  the  following,  symbols  and  formulas  used  in  photoelastic  analysis 
are  listed. 


Difference  in  principal  stresses; 


Kn 

t 


(Basic  photoelastic  equation) 
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where 


K 

- 

stress-optic  coefficient  in  lb/in. /order 

n 

XX 

fringe  order 

t 

a 

length  of  passage  of  the  light  vector 

0 *i  *  O2 

- 

stresses  in  1  and  2  directions 

For  the 

tension  specimen,  let  0^  “  0 

- 

P 

tv  (loaded  direction) 

Du 

where 

P 

- 

applied  load 

b 

m 

width  of  specimen  at  test  section 

h 

- 

thickness  of  specimen  ; 

then, 

K 

as 

n 

- 

2h 

P/bh  —  ,  t  -  2h  (for  reflection  polarscope) 
n 

P/n 

s 

slope  of  load  vs  fringe  order  curve 

s 

18.7  lb/in. /order  (as  determined  for  this  case) 

K 

m 

2P  2  x  18.7 

bn  0.5 

K 

= 

74.8  lb/ in. /order  for  this  type  of  resin 

With  the  stress-optic  coefficient  determined,  the  10-nil  glass  fiber 
slice  was  mounted  so  as  to  be  viewed  by  the  Budd  instrument. 

This  picture  under  polarized  light  is  shown  in  Figure  24.  Small 
grease  points  were  used  to  locate  spots  where  the  stress  measurements  were 
des ired . 

A  plot  of  the  difference  in  principal  stresses  at  a  radius  of  0.10  inch 
from  the  fiber  center  along  the  fiber  length  is  shown  in  Figure  25.  The 
0.10-inch  radius  was  chosen  for  plotting,  since  this  coincided  with  a  sta¬ 
tion  used  in  the  theoretical  analysis. 
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Loaded 


Figure  24.  Slice  of  10-mil  Glass  Fiber  in  Epoxy 
Resin.  Specimen  under  Reflection 
Polarscope  (Quarter  wave  plate  in 
system) 
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The  standard  Tardy  method  was  used  to  determine  the  fringe  order  at 
each  point  of  interest;  i.e.,  the  polarizer  and  analyzer  sections  of  the 
polarscope  were  oriented  each  time  to  bring  an  isoclinic  or  principal 
stress  direction  on  each  point.  The  analyzer  then  was  adjusted  to  bring 
the  lower  observed  fringe  order  to  the  point  of  interest.  The  fraction 
of  this  degree  movement  to  180°  gave  the  fractional  order  to  be  added  to 
the  lower  fringe  for  the  final  value. 

The  versatility  of  the  Budd  instrument  allows  an  additional  feature; 
i.e.,  to  pass  light  obliquely  into  the  specimen  and  thus  to  separate  the 
magnitude  of  the  stresses.  This  was  done  for  the  center  fiber  point  and 
only  for  the  specimen  without  external  load.  The  principle  of  oblique 
incidence  and  separation  of  the  principal  stresses  is  discussed  thoroughly 
in  Reference  1.  The  plot  of  magnitudes  of  these  stresses  for  the  center 
point  is  shown  in  Figures  26  and  27. 

The  theoretical  analysis  has  been  extended  to  the  case  where  the 
resin  surrounding  the  fiber  was  subjected  to  a  uniform  compressive  load 
in  the  direction  of  the  length  of  the  fiber. 

It  was  desired  to  accomplish  this  same  effect  in  the  experimental 

work. 


The  loading  scheme  for  this  operation  is  shown  in  Figure  28. 

A  minute  portion  at  both  ends  of  the  fiber  was  removed  inside  the 
resin  so  that,  in  effect,  it  was  not  loaded  under  a  pressure  application 
to  the  ends  of  the  specimen.  The  load  measurement  was  determined  by  the 
spring  load  cell  attachment  shown  in  Figure  28. 

The  results  of  this  determination  of  principal  stress  difference 
under  compressive  load  are  shown  also  in  Figure  25,  and  allow  a  direct 
comparison  of  the  unloaded  and  loaded  specimens . 

Unfortunately,  the  mirror  attachment  for  oblique  incidence  measure¬ 
ments  could  not  be  brought  close  enough  to  the  specimen  due  to  the  inter 
ference  of  the  load  cell.  Therefore,  it  was  not  possible  to  separate 
the  stresses  at  the  center  point. 
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Figure  27.  Axial  Tension  Stress  | °33 J  for  Center  Point  of 

the  Length  of  Fiber  versus  Radius  (Without 
External  Load) 


Figure  28.  Compressive  Loading  System  for 
the  Single  Glass  Fiber  Specimen 
(Fiber  diameter  -  0.010  inch) 
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MULTIFIBER  MODEL 

Glass  rods  0.275  inch  In  diameter  were  assembled  in  a  simulated 
cylinder  consisting  of  a  19-fibcr  array  and  an  epoxy  resin  of  the  same 
cycle  which  was  used  for  the  single-fiber  specimen. 

Figure  29  depicts  the  glass  rods  assembled  just  prior  to  the  surround 
ing  of  the  rods  with  resin. 

Upon  completion  of  the  fabrication  of  the  cylinder,  a  1/8- inch -thick 
slice  was  cut  from  this  specimen.  This  is  shown  in  Figure  30  under  polar¬ 
ized  light  (with  quarter  wave  plate  in  the  system). 

It  can  be  seen  that  the  specimen  shows  the  stresses  in  the  resin 
between  the  rods  to  be  symmetrical  in  nature. 

The  resin  between  a  group  of  three  rods  within  this  specimen  was 
selected  for  analysis  (the  selection  was  based  upon  uniformity  of  the 
distance  between  rods). 

Figure  31  shows  the  magnitude  of  the  difference  in  principal  stresses 
along  a  line  joining  the  center  of  the  resin  triangle  to  the  interface  of 
one  of  the  glass  rods. 

A  good  comparison  of  these  stress  magnitudes  was  obtained  with  the 
analysis  reported  in  Reference  2. 
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Figure  29.  Assembly  of  Multifiber  Glass 
Rod  Specimen  prior  to  Filling 
with  Resin 
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Figure  30.  1/8-inch-Thick  Slice  of  19-Rod 
Specimen  as  Viewed  under  Direct 
Polarscope  Optic  System 
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Figure  31.  Difference  In  Principal  Stresses  Along  Line 
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